Charge and Orbital Ordering and Spin State Transition 
Driven by Structural Distortion in YBaCo 2 Os 
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We have investigated electronic structures of antiferromagnetic YBaCo2 0s using the local spin- 
density approximation (LSDA) + U method. The charge and orbital ordered insulating ground 
state is correctly obtained with the strong on-site Coulomb interaction. Co 2+ and Co 3+ ions are 
found to be in the high spin (HS) and intermediate spin (IS) state, respectively. It is considered 
that the tetragonal to orthorhombic structural transition is responsible for the ordering phenomena 
and the spin states of Co ions. The large contribution of the orbital moment to the total magnetic 
moment indicates that the spin-orbit coupling is also important in YBaCo2 0s. 



PACS number: 75.50.Ee, 71.27. +a, 71.20.Be, 71.70.Ej 



Recently, an interesting spin state transition of the 
Co 2+ ion in YBaCo205 has been reported by Vogt et 
al. jjj using the neutron powder diffraction (NPD) mea- 
surements. The transition is induced by the long-range 
orbital and charge ordering of Co 2+ /Co 3+ ions. The or- 
dered oxygen-deficient double perovskite -RBaC^Os+a 
(R = rare-earths) has attracted much attention as a 
new spin-charge-orbital coupled system like manganites 
and also as a new Co-based colossal magnetoresistance 
(CMR) material. Indeed, giant magnetoresistance are 
observed for R=Gd and Eu, (po — Ph=7t)/po = 41% and 
40% for GdBaCo205 4 and EuBaCo205 4, respectively 

In the paramagnetic phase, YBaC^Os is crystal- 
lized in the tetragonal structure of the space group 
PA/mmm. It consists of double C0O5 square base pyra- 
midal backbone layers along the c-axis in which Y and 
Ba layers intervene alternatively and oxygens are defi- 
cient exclusively from the Y layers |^]-[|]. According to 
the valency consideration, Co 2+ and Co 3+ ions coex- 
ist similarly as the Mn 3+ /Mn 4+ covalency in hole-doped 
Lai-zSr^MnOs. Below T N ~ 330 K, YBaCo 2 5 under- 
goes a G-type antiferromagnetic (AFM) transition and 
the lattice changes slightly from the tetragonal to or- 
thorhombic structure At Tqo ~ 220 K, a pronounced 
upturn is observed in the resistivity indicating that an- 
other transition takes place, i.e., the long-range charge 
and orbital ordering of Co 2+ /Co 3+ ions. The stripe type 
charge ordering is formed in the ab plane. Further, the 
spin state of the Co 2+ ion changes from the low to high 
spin across Tqo, which is evidenced by the increased mag- 
netic moment per Co ion from 2.10/zb at 300 K to 3.45/^b 
at 25 K. More recently, for the isostructural HoBaCo20s, 
essentially the same features are observed of Tn ~ 340 K 
and T co ~ 210 K ||. 

The phenomenon of the spin state transition is 
observed usually in cobaltates such as LaCo03 and 
Lai-ajSr^CoOs. In LaCoOs, the magnetic ground state 
of Co 3+ ion corresponds to the low spin (LS) state with 
t% g e° g . Upon heating, a successive spin state transition to 
an intermediate spin (IS) state and then to a high spin 
(HS) state occurs f|. Note that, distinctly from the case 



in LaCoOs, the spin state transition in YBaCo205 seems 
to occur for Co 2+ ion. In fact, this issue is still contro- 
versial. Based on the reduction of the magnetic suscepti- 
bility below 220K, Akahoshi and Ueda suggested that 
the AFM transition takes place at T ~ 220 K which is 
associated with a spin state transition of Co 3+ from the 
HS to LS state upon cooling. Thus, the nature of the 
spin state transition and the interplay of the spin state 
with the charge and orbital ordering are still unclear. 

To reveal the mechanism of the spin state transition as 
well as the charge and orbital ordering, we have explored 
the electronic structure of the G-type AFM YBaCo2C-5 
using the local-spin density approximation (LSDA) + U 
scheme implemented in the linearized muffin-tin orbital 
band method §,§. 

We have employed two structural data of nearly tetrag- 
onal structure at 300 K (LI) and orthorhombic structure 
at 25 K (L2) Q]. In the LI structure, all the Co sites 
have an equal average bond length of d(Co-O) = 1.97 A. 
Whereas in the L2 structure, there are two different kinds 
of Co sites, Col and Coll. The bond length at Col sites 
is extended to d(CoI-O) = 2.03 A and that at Coll sites 
becomes shortened to d(CoII-O) = 1.92 A. The G-type 
AFM spin order is assumed in all our calculations. 

In Fig. [|, we have compared the LSDA Co 3c? par- 
tial density of states (PDOS) for the LI and L2 struc- 
tures. Most notable is the change of the exchange split- 
ting which becomes larger for Col site and smaller for 
Coll in L2 than that for Co in LI. This indicates that, 
due to the lattice distortion, 3d electrons in L2 becomes 
more localized for Col and less localized for Coll in com- 
parison to those for Co in LI. The calculated spin mag- 
netic moments are ps = 1.68/is and 0.95/is for each Col 
and Coll site, respectively, in L2 and ps = \.1?>pb for 
Co in LI. The sizes of the spin moments are consistent 
with the degree of 3d electron localization. Hence, the 
LSDA gives a qualitatively good information about the 
structural transition effects on Co 3d electron states. The 
LSDA, however, yields an incorrect metallic ground state 
which pertains even after the structural transition from 
the LI to L2 structure. This is different from the exper- 
iment which shows unambiguously the semiconducting 
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FIG. 1. The LSDA PDOS of Co 3d electrons in the LI 
(300 K phase) and L2 (25 K phase) structures. In L2, there 
are two different types of Co sites, Col and Coll. Neighboring 
Col and Coll sites are polarized antiferromagnetically. As 
manifested by the exchange splitting, 3d electrons are more 
localized for Col site and less localized for Coll than those for 
Co in LI. 

resistivity behavior below 220 K H,|j7). 

Using the LSDA, one cannot expect proper description 
of localized Co 3d electrons in YBaC^Os. To resolve the 
above problem, we have applied the LSDA + U method 
with parameter values of U = 5.0 eV and J = 0.89 eV. 
Although there is an arbitrariness of the J7-value in our 
calculation, the LSDA + U results are usually not much 
sensitive to the used U- value within AU ~ ±1 eV |Lo[ . 
The spin-orbit interaction is also taken into account in 
the self-consistent variational loop, because Co 3c? elec- 
trons are expected to retain atomic properties to some 
extent due to their localized nature. 

Figures | and | present Col and Coll 3d PDOS, re- 
spectively, in the L2 structure obtained by the LSDA + 
U calculations. In the bottom panels, t2 g and e g decom- 
positions of 3d PDOS are also provided. It is amusing 
to note that the energy gap of E g ~ 0.6 eV opens at the 
Fermi level E-p and so YBaCo20s becomes an insulator 
as expected in consideration of the Co 2+ /Co 3+ charge 
ordering. For Col 3d electrons, the majority-spin bands 
are fully occupied by three electrons in ti g and two elec- 
trons in e g bands. The minority-spin t2 g bands are only 



partially occupied by two electrons, and one ti g and two 
e g bands are almost empty. Considering the pyramidal 
environment of Co, the occupied t2 g states corresponds 
to d zx and d yz while the empty t2 g to d xy . Accordingly, 
the nominal valency and the 3d electron configuration 
at Col site are assigned to be Co 2+ and 3d 7 (t^e 2 ,), re- 
spectively pl[ . Hence, Co 2+ ion is in the HS state with 
spin magnetic moment of us = 3[Ib {S = 3/2), which is 
consistent with the NPD data. 

For Coll 3d electrons, the majority-spin bands are not 
fully occupied with a split-off empty e g state above E-p 
(see Fig. ||). For the minority-spin bands, the situation 
is similar to the case of Col. It is thus possible to identify 
the 3d electron configuration of Coll as 3d 6 (t| g e g) an d 
the valency as Co 3+ . With one less electron than Co 2+ , 
the lower d 3z 2_ r 2 out of two e g states in the majority-spin 
bands is occupied and the upper d x i_ y i becomes empty. 
Hence, the spin state of Co 3+ ion is the IS state which 
is also in agreement with the experimental analysis of 
fi s = 2/is (S = 1). 

The calculated charge occupancies of each Co 3d or- 
bitals are shown in Table |. For Co 2+ ion, the majority- 
spin bands are almost completely occupied by t2 g ]= 
2.98 and e g f= 1.99 electrons, while the minority- 
bands are only partially occupied by t2 g 1= 2.02 and 
c g 1= 0.37 electrons. For Co 3+ ion, it is noticeable 
that the majority-spin e g states are partially occupied 
by e g I— 1.54 electrons. Hence, the calculated charge 
occupancies are consistent with the nominal valencies of 
Co 2+ and Co 3+ , if one takes into account the band hy- 
bridization effects. 

Vogt et al. JIJ have deduced magnetic moments from 
the NPD experiments as /Lt cxp = 4.2/zb and 2.7/ib for 
each Co 2+ and Co 3+ ion, respectively. In the analy- 
sis of the experiments, they counted only the spin mo- 
ment contribution, assuming that the orbital moment is 
quenched. However, the orbital moment is only partially 
quenched in YBaCo20s. In Table |, we have summarized 
the calculated magnetic moments using the LSDA + U 

TABLE I. The calculated charge occupancies of Co 3d or- 
bitals (electrons) and magnetic moments (/is) of YBaCo2 0s 
in the LSDA + U method. 
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6 HoBaCo205 in Ref. ||. The assignment of the mag- 
netic moments to each Co ion is corrected (see text). 
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2. The LSDA + U PDOS of Col 3d electrons in the 



L2 (25 K phase) structure. In the bottom panel, ti a and e 9 
decompositions of 3d PDOS are also provided. The major- 
ity-spin bands are fully occupied, while in the minority-spin 
bands, only 2/3 states of t2 g bands are occupied. The nomi- 
nal valency of Col site is Co 2+ with 3d 7 (t 
in the HS state. 
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FIG. 3. The LSDA + U PDOS of Coll 3d electrons in 
the L2 (25 K phase) structure which is antiferromagnetically 
polarized to that of neighboring Col. Unoccupied e g states 
in the majority-spin bands are clearly visible above Ef. The 
valency and the 3d electron configuration of Coll site is nom- 
inally identified as Co 3+ and 3d 6 (t 
ion is in the IS state. 



el), respectively. Co 3+ 



method. For Co 2+ ion, the spin and orbital magnetic 
moments are /is = 2.61/is and = 1.04/is, respec- 
tively, and for Co 3+ ion, fi$ = 1.84/is and fj,L = 0.40/^s. 
The orbital moment of Co 2+ ion is as much as that of 
CoO jlj,[ll|, and the orbital moment of Co 3+ is com- 
parable to that of NiO 0@. The non- negligible or- 
bital moment, which originates from the localized na- 
ture of Co 3d electrons, suggests that YBaC^Os should 
fall in a class of strongly correlated electron system like 
CoO and NiO. The calculated total magnetic moments 
of jtxtot = and 2.24^b for Co 2+ and Co 3+ ion, 

respectively, are only slightly smaller than the experi- 
mental values. Evidently, this interpretation will also be 
valid for HoBaCo20s. Suard et al. || have improperly 
assigned the NPD measured /i oxp = and 2.7/jb in 

HoBaCo 2 05 to Co 3+ and Co 2+ ion, respectively, How- 
ever, their assumptions of spin-only moments and the HS 
states for both Co 2+ and Co 3+ ions are discarded by the 
present results. To test the calculated results, more ex- 
perimental works like the x-ray scattering measurement 
are encouraged, in which separate determination of the 
spin and orbital moments are possible. 



In Fig. |4|, we have plotted the geometry of the orbital 
ordering which is obtained from the orbital dependent 
occupancy of the 3d minority-spin states at each Co site. 
At Co 2+ sites, the orbitals are aligned along a-axis, while 
at Co 3+ sites, the orbitals are along o-axis. This feature 
is understandable by considering that the bond length of 
<i(Co 2+ -0) is larger in a-axis than in o-axis and vice versa 
for that of d(Co 3+ -0). As for the charge ordering configu- 
ration, Co 2+ and Co 3+ chains of a stripe type are formed 
in the ab plane along o-axis, which are alternating in the 
a and c direction [QJ|]. This is in contrast to the charge 
ordering observed in the isostructural YBaM^Os. In 
YBaMn 2 5 , the Mn 2+ /Mn 3+ orders in a checkerboard 
type |H| . This difference gives rise to the different mag- 
netic structures: G-type AFM phase for YBaCo20s and 
G-type ferrimagnetic phase for YBaM^Os. The theo- 
retical result of the charge and orbital ordering geome- 
try coincides with the experimentally proposed one jl],|| . 
Thus, it can be inferred that the deformed bond lengths 
of d(Co-O) determine the charge and orbital ordering ge- 
ometry. 

The G-type AFM ordering in YBaCo20s is consistent 
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FIG. 4. The orbital ordering geometry of the occupied 
minority-spin states at each Co site. At Co 2+ sites, 3d orbitals 
are aligned along the a-axis while at Co 3+ sites, along the 
b- axis. 



with the above charge and orbital ordering geometry. In 
the Co 3+ chains, the kinetic-exchange energy gain be- 
tween the occupied d yz and d zx states and the empty 
d x 2_ y 2 state for the AFM configuration of neighboring 
sites stabilizes the AFM spin ordering. In a similar way, 
the AFM ordering between neighboring Co 3+ and Co 2+ 
ions can be explained. The AFM ordering in the Co 2+ 
chains, however, is hard to understand in terms of the di- 
rect kinetic-exchange gain, because the overlap integral 
between two neighboring Co 2+ ions would be negligible 
as seen in Fig. [I| Instead, the AFM interaction in the 
Co 2+ chains is expected to be derived indirectly via the 
Co 2+ -Co 3 + and Co 3+ -Co 3+ AFM interactions. 

As mentioned above, the structural transition plays a 
crucial role in determining the ground state properties 
of YBaCo205. Although the tetragonal to orthorhom- 
bic structural transition occurs simultaneously with the 
G-type AFM transition at Tn ~ 330 K, the structural 
deformation is not significant above 220 K. Only near 
Tqo ~ 220 K, the lattice splitting between a- and &-axis 
becomes pronounced and the charge and orbital ordering 
emerges with the Co 2+ spin state transition from the low 
to high spin. Furthermore, it is known that the long- 
range charge ordering and the spin state transition are 
very sensitive to the oxygen stoichiometry There- 
fore, the structural distortion is thought to be respon- 
sible for the orderings and the spin state transition by 
inducing different local environment at each Co ion site. 
This feature implies that the electron-lattice interaction 
is very important in this system. A detailed study on 



the electron-phonon interaction effects in YBaC^Os is 
urgently demanded. 

In conclusion, we have performed the LSDA + U calcu- 
lations for a new spin-charge-orbital-lattice coupled sys- 
tem YBaCo 2 05. It is found that the Co 2+ /Co 3+ charge 
and orbital ordering and the Co 2+ HS state transition 
are closely correlated with the lattice distortion from the 
tetragonal to orthorhombic structure. The orbital mo- 
ment has a substantially large contribution to the total 
magnetic moment. All of the effects of the Coulomb cor- 
relation, the spin-orbit coupling, and the electron-phonon 
interaction should be properly taken into account to un- 
derstand physical properties of YBaCo205. 

Acknowledgements— This work was supported by the 
KOSEF (1999-2-114-002-5) and by the Brain Korea 21 
Project. 



[1 

[2 
[3. 
[4 
[5 
[6 

[r 

[8 
[9 

[io; 
in 

[12 
[13 
[14 

[is; 



T. Vogt, P. M. Woodward, P. Karen, B. A. Hunter, P. 
Henning, and A. R. Moodenbaugh, Phys. Rev. Lett. 84, 
2969 (2000). 

W. Zhou, C. T. Lin, and W. Y. Liang, Adv. Mater. 5, 
735 (1993); W. Zhou, Chem. Mater. 6, 441 (1994). 

C. Martin, A. Maignan, D. Pelloquin, N. Nguyen, and B. 
Raveau, Appl. Phys. Lett. 71, 1421 (1997). 

A. Maignan, C. Martin, D. Pelloquin, N. Nguyen, and B. 
Raveau, J. Solid State Chem. 142, 247 (1999). 

E. Suard, F. Fauth, V. Caignaert, I. Mirebeau, and G. 
Baldinozzi, Phys. Rev. B 61, Rll 871 (2000). 

M. A. Korotin, S. Yu. Ezhov, I. V. Solovyev, V. I. Anisi- 
mov, D. I. Khomskii and G. A. Sawatzky, Phys. Rev. 54, 
5309 (1996); and references therein. 

D. Akahoshi and Y. Ueda, J. Phys. Soc. Jpn. 68, 736 

(1999) . 

S. K. Kwon and B. I. Min, Phys. Rev. Lett. 84, 3970 

(2000) . 

A. I. Liechtenstein, V. I. Anisimov, and J. Zaanen, Phys. 
Rev. B 52, R5467 (1995). For review, see V. I. Anisi- 
mov, F. Aryasetiawan, and A. I. Liechtenstein, J. Phys. 
: Condens. Matter 9, 767 (1997). 

S. K. Kwon and B. I. Min, " Unquenched large orbital 
magnetic moment in NiO", Phys. Rev. B (in press). 
The valency assignments are nominal, because the num- 
ber of occupied d-electrons inside the atomic sphere of 
Co is not an integer due to the hybridization with other 
bands (see Table I). 

A. Svane and O. Gunnarsson, Phys. Rev. Lett. 65, 1148 
(1990). 

I. V. Solovyev, A. I. Liechtenstein, and K. Terakura, 
Phys. Rev. Lett. 80, 5758 (1998). 

V. Fernandez, C. Vettier, F. de Bergevin, C. Giles, and 
W. Neubeck, Phys. Rev. B 57, 7870 (1998). 

F. Millange, E. Suard, V. Caignaert, and B. Raveau, 
Mater. Res. Bull. 34, 1 (1999). 



4 



